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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development  Center 
(AEDC),  Air  Force  Systems  Command  (AFSC),  at  the  request  of  the  Arnold  Engineering 
Development  Center  (AEDC/DOT).  The  AEDC/DOT  project  manager  was  Dr.  Keith 
Kushman.  The  results  were  obtained  by  Calspan  Corporation/ AEDC  Division,  operating 
contractor  for  the  aerospace  flight  dynamics  testing  effort  at  the  AEDC,  AFSC,  Arnold  Air 
Force  Station,  Tennessee.  Analysis  of  the  results  was  completed  December  1,  1985.  The  effort 
was  conducted  under  Air  Force  Project  Number  DA04PW.  The  manuscript  was  submitted 
for  publication  on  March  28,  1986. 
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1.0  INTRODUCTION 

The  nozzle  afterbody  flow  field  is  a  complex  interaction  between  a  jet  exhaust  plume, 
the  local  free-stream  flow,  and  the  model  afterbody  surfaces.  Variables  that  affect  the  afterbody 
flow  field  are  nozzle  geometry,  afterbody  geometry,  nozzle  afterbody  base,  exhaust  jet 
temperature,  exhaust  jet  gas  constant,  exhaust  jet  specific  heat  ratio,  and  nozzle  pressure 
ratio  (NPR).  Other  considerations  are  model  scale,  support  system  interference,  downstream 
hardware  interference,  and  free-stream  effects  such  as  Mach  number,  Reynolds  number,  and 
humidity.  Much  work  has  been  performed  at  AEDC  directed  toward  understanding  these 
effects  and  manipulating  them  in  a  controlled  environment  to  correctly  simulate  flight 
conditions  in  the  wind  tunnel.  Kennedy  in  Ref.  1  summarized  AEDC’s  work  in  this  area 
up  to  the  year  1980.  Later  investigations  are  reported  by  Price  and  Lucas  in  Refs.  2  and 
3,  respectively.  With  the  increase  in  capabilities  of  mainframe  computers,  computational  fluid 
dynamics  (CFD)  presents  the  opportunity  to  develop  a  tool  for  flow  analysis  complementary 
to  the  wind  tunnel.  Accurate  flow  simulation  by  CFD  methods  can  provide  details  that  are 
difficult,  if  not  impossible,  to  obtain  in  wind  tunnels.  Information  such  as  flow  field  velocity 
vectors  and  pressure  contours  in  addition  to  typical  wind  tunnel  data  such  as  model  surface 
pressure  distributions  and  forces  can  be  obtained  from  numerical  solutions.  Long-term 
objectives  are  for  CFD  to  vastly  decrease  the  requirements  for  nozzle  afterbody  wind  tunnel 
testing  through  numerical  simulation;  however,  presently,  computer  size  and  computational 
cost  are  limiting  factors. 

It  was  demonstrated  by  Jacocks,  Peters,  and  Guyton  in  Ref.  4  that  CFD  can  be  a  valuable 
tool  for  use  on  the  nozzle  afterbody  problem.  In  this  work,  an  axisymmetric,  Navier-Stokes 
code  (Ref.  5)  was  used  to  calculate  the  flow  about  the  AGARD  15-deg  boattail  afterbody 
with  a  full  flowing  jet.  The  solutions  generated  were  compared  to  wind  tunnel  test  data.  Figures 
1  and  2  show  details  of  the  model  used  to  obtain  the  experimental  data  (Ref.  6)  and  Fig. 

'  3  shows  a  single  plane  view  of  the  axisymmetric  computational  mesh.  Unlike  the  computational 
model,  the  experimental  article  had  a  finite  nozzle  base  thickness  for  instrumentation  purposes. 
Shown  in  Fig.  4  is  a  comparison  of  the  measured  and  computed  pressure  distributions  over 
the  afterbody  at  Mach  0.9  and  design  NPR.  The  two  results  agree  well,  except  in  the 
recompression  region  where  the  predicted  results  indicate  a  higher  presure  than  was  measured 
at  the  end  of  the  afterbody.  This  difference  may  be  attributed  to  the  fact  that  there  was  no 
model  base  present  in  the  numerical  solution.  Figure  5  presents  a  comparison  of  experimental 
and  computational  pressure  drag  versus  NPR.  The  shape  and  magnitude  of  the  variations 
are  similar.  The  computations  correctly  predict  that  with  pressure  ratio  increasing  from  the 
jet-off  condition,  an  initial  decrease  in  drag  occurs,  followed  by  an  increase  in  drag  to  a 
maximum  jet-on  ’‘peak  drag”  value.  This  region  is  characterized  by  an  overexpanded  jet 
plume  and  is  referred  to  as  the  “drag  bucket."  Further  increases  in  NPR  cause  a  continuing 
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Tunnel 


Figure  1.  Axisymmetric  model  dimensions  and  location  in  tunnel  16T. 
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Figure  2,  Afterbody  contour  and  surface  coordinates  for  axisymmetric  model. 
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Figure  3.  Nozzle  afterbody  computational  mesh  for  axisymmetrfc  Navier-Stokes  calculations. 

decrease  in  drag  produced  by  the  increasing  blockage  of  the  underexpanded  jet  plume.  The 
computed  values  of  drag  coefficient  are  higher  than  experiment  in  the  overexpanded  jet  regime; 
however,  computational  results  correctly  predict  drag  coefficient  values  at  underexpanded 
jet  conditions.  The  disagreement  in  the  “drag  bucket”  region  may  also  be  attributable  to 
no  base  being  present  in  the  computations,  causing  a  larger  influence  of  jet  entrainment  on 
the  afterbody  that  was  actually  experienced.  Figure  6  presents  Mach  number  contour  plots 
from  the  Navier-Stokes  computations.  The  effects  of  increasing  pressure  ratio  are  easily 
observed  in  the  plume  Mach  number  contours;  however,  the  effects  causing  variation  in  the 
integrated  drag  presented  in  Fig.  5  are  not  as  evident  when  the  afterbody  region  is  viewed. 

The  effects  of  nozzle  geometry,  jet  temperature,  and  Mach  number  were  also  investigated 
computationally  by  Jacocks,  et  al.  The  axisymmetric,  Navier-Stokes  code  provided  accuracy 
comparable  to  the  previous  nozzle  pressure  ratio  analysis  for  these  variables  as  well.  In  general. 
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Figure  4.  Computational  and  experimental  pressure  distributions  for  axisymmetric 
Navier-Stokes  calculations. 

this  investigation  demonstrated  the  utility  and  capability  of  a  Navier-Stokes  code  for  the  study 
of  the  nozzle  afterbody  problem. 

The  high  cost  of  computer  resources  required  for  CFD  solutions  makes  it  desirable  to 
investigate  the  possibilities  of  using  the  most  simplified  equations  available  that  will  correctly 
predict  the  flow  field  of  the  particular  problem  under  study.  The  purpose  of  this  project 
was  to  determine  the  feasibility  of  using  a  viscid/inviscid  interaction  code  to  compute  three- 
dimensional  nozzle  afterbody  flow  fields.  It  was  believed  that  some  of  the  shortcomings  of 
previous  attempts  at  solving  this  problem  might  be  attributed  to  use  of  potential  flow  for 
the  inviscid  calculations,  solid  plume  simulators  that  do  not  account  for  jet  mixing,  or  direct 
integral  boundary  layers  that  do  not  allow  flow  separation.  The  code  developed  during  this 
investigation,  hereafter  referred  to  as  AR02P,  solves  the  Euler  equations  for  nozzle  afterbody 
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=  0.9 

Re  =  1.5  x  106 


Figure  5.  Effect  of  nozzle  pressure  ratio  on  integrated  drag  for  axi  symmetric  Navier-Stokes 
calculations. 

flow  fields  and  the  real  jet  plume  and  uses  an  inverse  boundary-layer  method  for  the  viscid 
correction. 

A  nominal  value  of  the  measurement  uncertainty  of  nozzle  afterbody  data  is  about  five 
aircraft  drag  counts  (Ref.  7).  Experience  has  indicated  that  one  can  expect  better  repeatability 
than  five  counts  for  data  obtained  by  pressure  integration  and  possibly  greater  than  five  counts 
for  data  obtained  with  a  force  and  moment  balance.  Five  aircraft  drag  counts  on  a  typical 
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fighter  aircraft  is  comparable  to  approximately  75  drag  counts  based  on  a  body  cross-sectional 
area  such  as  was  used  in  Ref.  5.  If  one  is  attempting  to  complement  tunnel  data  by  filling 
in  points  in  a  matrix  or  by  extrapolating  data  by  CFD  analysis,  the  accuracy  of  the  numerical 
solutions  should  be  at  least  as  good  as  the  tunnel  data.  Therefore,  the  success  criteria  for 
evaluating  the  code  capabilities  are  that  calculations  be  within  75  body  drag  counts  of  the 
experiment.  This  is  not  to  say  the  absolute  value  of  the  CFD  solution  must  be  within  75  body 
drag  counts,  but  rather,  incrementally  the  data  should  fall  within  a  150  drag  count  band 
about  a  line  of  the  same  slope  as  the  data.  If  one  is  investigating  an  unknown  situation,  for 
example  correcting  tunnel  data  for  some  effect,  then  any  computations  that  give  solutions 
closer  to  the  real  answers  are  beneficial. 

To  determine  the  utility  of  the  CFD  code  for  nozzle  afterbody  analysis,  a  check  case  with 
known  answers  was  used.  For  the  aforementioned  Navier-Stokes  code  and  for  the  ensuing 
computational  results,  the  AGARD  15-deg  boattail  afterbody  provides  a  simple,  but 
representative  model  with  a  sufficient  wind  tunnel  database.  Experimental  data  obtained  on 
the  AGARD  afterbody  were  on  a  strut-mounted  model  (Fig.  1).  The  strut  was  not  modeled 
for  the  Navier-Stokes  calculations,  nor  was  it  included  in  the  present  computational  grid. 
Interference  from  the  strut  no  doubt  had  some  influence  on  the  experimental  data,  although 
the  magnitude  of  this  effect  has  not  been  computed.  At  a  given  Mach  number,  the  effect 
of  the  strut  interference  on  the  drag  variation  with  NPR  is  believed  to  be  very  small.  Therefore, 
comparing  slopes  of  the  drag  variation  with  NPR  between  the  experiment  and  the  CFD  code 
is  believed  to  be  valid.  Should  the  code  prove  capable  of  computing  the  effects  of  parameters 
varied  in  the  experiments  such  as  NPR,  Mach  number,  or  jet  temperature  for  the  check  case, 
then  it  is  reasonable  to  expect  it  to  do  the  same  for  more  complicated  geometries.  However, 
using  the  code  to  predict  effects  not  verified  by  experiment  would  be  subject  to  question. 

2.0  CODE  DEVELOPMENT  AND  ANALYSIS 
2.1  AR02P  COMPONENTS 

The  inviscid  core  for  AR02P  is  the  AROl  code  developed  by  Jacocks  and  Kneile  (Ref. 
8).  It  was  created  to  solve  the  three-dimensional,  unsteady  Euler  equations  in  Cartesian 
coordinates  using  a  finite  volume  (volume  flux)  approach.  The  basic  numerical  algorithm 
is  the  explicit  predictor-corrector  scheme  of  MacCormack  (Ref,  9).  The  viscous  component 
is  a  boundary-layer  correction  routine  developed  by  Whitfield,  Swafford,  and' Jacocks  (Ref. 
10).  It  is  an  inverse  boundary-layer  method  that  calculates  turbulent  boundary  layers  with 
and  without  separation.  “Inverse”  refers  to  the  technique  whereby  the  boundary-layer 
displacement  thickness  is  prescribed,  as  opposed  to  the  direct  method  of  calculation  in  which 
the  pressure  distribution  is  prescribed.  The  effect  of  the  boundary  layer  is  represented  in 
the  inviscid  code  by  a  transpiration  condition  at  the  appropriate  boundaries. 


11 


AEDC-TR-86-18 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

X/L 

Figure  7.  Computational  and  experimental  pressure  distributions  for  axisymmetric  Euler, 
viscid/inviscid  calculations. 


Calculations  were  made  to  verify  the  capability  of  AR02P  prior  to  addition  of  the  exhaust 
jet.  Figure  7  shows  the  variation  in  pressure  distribution  over  an  F-15  equivalent  body  of 
revolution.  The  solutions  were  generated  with  the  code  operating  in  the  viscid/inviscid  mode. 
As  with  the  experimental  model,  a  solid  sting  was  included  in  the  calculations;  there  is  no 
jet  flow.  The  variation  of  pressure  coefficient  is  in  good  agreement  with  the  measured  values. 

The  code  was  completed  by  the  addition  of  a  jet  plume  calculation  to  the  inviscid  routine. 
The  jet  exhaust  model  is  based  on  conical  flow  and  allows  variation  in  NPR,  temperature 
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ratio,  nozzle  area  ratio,  and  nozzle  exit  angle.  The  addition  of  the  plume  caused  a  significant 
instability  in  the  flow  field  calculation  at  the  trailing  edge  of  the  afterbody  where  the  free- 
stream  and  jet  exhaust  flow  meet.  To  stabilize  the  calculations  a  new  smoothing  scheme 
developed  by  J.  L.  Jacocks  and  K.  C.  Reddy  was  integrated  into  the  program.  The  routine 
is  based  on  a  total  variation  diminishing  technique  originated  by  Stephen  Davis  (Ref.  1 1). 

2.2  COMPUTATIONAL  GRID 

The  computational  grids  used  in  CFD  are,  in  general,  problem  unique,  varying  with  model 
geometry  to  adequately  distribute  points  to  resolve  expected  gradients.  However,  one  can 
consistently  use  the  same  type  of  grid  for  the  same  type  of  problem.  In  fact,  it  is  desirable 
to  do  so  if  the  information  generated  by  the  different  calculations  will  be  used  together  for 
comparison  analysis. 

The  type  of  grid  chosen  for  AR02P  is  a  single  piece  mesh  with  a  sectioned  plume.  Figure 

8  depicts  the  transition  of  the  grid  from  the  computational  plane  to  the  physical  plane.  The 
top  drawing  shows  the  computational  plane  with  two  distinct  sections  outlined.  The  larger 
contains  the  free-stream  flow  and  the  smaller  section  the  jet  plume  flow.  The  transition  of 
the  free-stream  flow  regime  to  the  physical  plane  shows  little' change,  but  the  plume  section 
is  rotated  90-deg  clockwise,  moved  underneath  the  free-stream  grid,  and  overlapped  one  grid 
line.  Boundary  conditions  are  applied  over  the  overlapped  row  of  the  plume  grid  to  provide 
a  transparent  connection  between  the  plume  and  free-stream  grid  sections.  Also  shown  in 
Fig.  8  are  representations  of  the  boundary  conditions  in  the  physical  plane  associated  with 
each  of  the  mesh  boundaries.  The  top  boundary  is  specified  to  simulate  a  porous  wind  tunnel 
wall.  Other  conditions  represented  are:  free-stream  and  jet  exhaust  inflow,  downstream  out¬ 
flow,  a  solid  wall  for  the  model  surface  modified  by  boundary-layer  transpiration,  and  a 
reflection  plane  at  the  plume  centerline.  Inital  conditions  are  uniform  flow  at  a  specified 
free-stream  Mach  number  and  a  boundary-layer  displacement  thickness  distribution 
corresponding  to  a  cylindrical  afterbody  at  a  specified  Reynolds  number.  Calculations  were 
initiated  26.5  in.  from  the  base  of  the  model,  which  is  approximately  10  in.  upstream  of 
the  boattail  region.  Initial  boundary-layer  displacement  thicknesses  at  the  beginning  of  the 
calculation  were  0.19  in.  and  0.17  in.  for  Mach  numbers  0.6  and  0.9,  respectively.  Figure 

9  shows  the  physical  plane  mesh  used  for  the  AGARD  15-deg  boattail.  Notice  the  coarseness 
of  the  mesh  in  the  vertical  direction  compared  to  the  Navier-Stokes  code  mesh  presented 
in  Fig.  3.  Since  the  Euler  equations  do  not  have  the  viscosity  terms,  no  boundary  layer  is 
calculated.  Therefore,  there  is  no  requirement  for  a  dense  grid  at  the  surface  to  resolve  the 
boundary-layer  gradients.  The  boundary  layers  used  with  the  inviscid  codes  such  as  AR02P 
merely  provide  a  correction  to  the  boundary  conditions  specified  on  the  body  surface,  but 
no  boundary  layer  is  calculated  as  part  of  the  flow  field. 
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Computational  Plane 


Figure  8.  Transformation  of  the  Euler  code  grid  from  the  computational  plane  to  the  physical 
plane,  In  the  X-Y  plane. 
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2.3  CODE  VALIDATION 

Calculations  were  made  of  the  flow  field  about  the  AGARD  body  using  AR02P  with 
the  boundary-layer  calculation  disabled,  that  is,  inviscid.  Figure  10  shows  the  variation  in 
pressure  distribution  over  the  afterbody  at  the  design  NPR  condition  and  a  free-stream  Mach 
number  of  0.6.  The  computed  distribution  predicts  a  lower  minimum  and  a  higher  maximum 
pressure,  indicating  that  the  calculation  senses  the  actual  afterbody  shape  with  no  distortion 
from  viscous  effects.  Further  calculations  were  made  over  a  range  of  NPR’s.  A  comparison 
of  the  integrated  drag  from  the  solutions  with  experiment  is  presented  in  Fig.  II.  The  slope 
of  the  inviscid  calculations  is  in  the  right  direction  but  is  greater  than  that  of  the  measured 
data.  One  would  expect  a  difference  in  slopes  in  that  direction  considering  there  are  no  viscous 
effects  in  the  calculations.  Considering  the  previously  defined  accuracy  requirements,  the 
predictions  do  provide  information  concerning  the  direction  of  the  integrated  drag  variation 
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Figure  10.  Inviscid  pressure  distribution  from  Euler  code  at  design  NPR. 

with  NPR,  but  they  are  not  within  a  150  drag  count  band  for  the  NPR  range  under  study. 
Thus,  a  viscous  correction  is  required  to  make  AR02P  useful  for  nozzle  afterbody  analysis. 

There  are  many  variables  in  CFD  codes  that  affect  the  solution,  yet  have  no  physical 
meaning  with  respect  to  a  real  flow  field.  Three  such  variables  that  are  used  to  “fine  tune” 
the  AR02P  code  are:  (1)  the  Courant-Friedrich-Lewy  stability  criterion  (CFL  number)  that 
determines  the  speed  with  which  the  numerical  solution  progresses,  (2)  OMEGA  -  a  parameter 
that  affects  the  amount  of  smoothing  applied  to  the  flow  field,  and  (3)  NV  -  the  number 
of  inviscid  iterations  between  boundary-layer  corrections. 

The  CFL  number  can  vary  between  zero  and  one,  where  one  represents  the  maximum 
time  stepping  speed.  Several  solutions  were  calculated  over  the  CFL  range  for  Mach  0.6  and 
design  NPR  to  determine  the  best  stability  condition.  Figure  12  shows  the  variation  of  pressure 
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0.04 


Figure  11.  Experimental  and  computational  variation  of  NPR  with  afterbody  drag  for  the 
inviscid  Enler  code. 

distribution  for  two  of  the  solutions,  CFL  =  0.6  and  CFL  =  0.9.  CFL  =  0.6  produced 
the  most  stable  solution  and  was  thus  used  throughout  the  rest  of  the  investigation. 

A  value  of  zero  for  the  smoothing  parameter  OMEGA  represents  the  maximum  amount 
of  smoothing  available,  while  any  number  greater  than  zero  is  a  decrease  in  smoothing.  It 
is  desirable  to  use  as  little  smoothing  as  possible  in  that  smoothing  introduces  error  into  the 
solution.  Figure  13  presents  the  variation  in  pressure  distribution  for  two  values  of  OMEGA, 
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Figure  12.  Variation  of  pressure  distribution  with  the  CFL  stability  criterion. 

zero  and  one.  {For  reference,  an  OMEGA  as  high  as  four  is  not  uncommon  for  other 
problems.)  It  is  obvious  that  OMEGA  =  0  produces  a  solution  closer  to  the  experimental 
data.  The  result  was  not  unexpected  considering  the  instabilities  encountered  earlier  with 
the  trailing  edge.  Consequently,  the  maximum  smoothing  available,  OMEGA  =  0,  was  used 
for  all  subsequent  calculations. 

Correcting  the  inviscid  solution  for  viscous  effects  on  each  iteration  would  be  very  time 
consuming.  Fortunately,  it  is  not  required.  Typically,  the  inviscid  part  of  a  code  such  as  AR02P 
is  run  for  a  set  number  of  iterations  (NV)  and  then  receives  information  for  a  viscous  correction 
to  its  boundary  conditions  from  a  boundary-layer  routine.  The  code  is  run  inviscidly  for 
NV  time  steps  again,  and  another  viscous  correction  is  made  based  on  the  latest  inviscid  solution 
and  so  forth.  Solutions  were  calculated  with  AR02P  at  Mach  0.6  and  design  NPR  for  various 
NV’s.  Figure  14  shows  solutions  for  NV  =  40  and  NV  =  80.  The  difference  between  the 
two  pressure  distributions  is  minimal,  although  NV  =  40  does  appear  to  be  slightly  more 
stable.  Therefore,  NV  =  40  was  chosen  for  the  rest  of  this  analysis. 
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Figure  13.  Variation  of  pressure  distribution  with  smoothing  parameter,  G. 

Figure  IS  presents  the  data  shown  in  Fig.  13  with  the  viscid/inviscid  solutions.  Again, 
the  calculated  values  show  a  decrease  in  CD  with  increasing  NPR.  However,  rather  than 
duplicating  the  slope  of  the  experimental  data,  the  slope  of  the  viscid/inviscid  calculations 
unexpectedly  is  steeper  than  the  inviscid  calculations.  Three  areas  were  identified  that  might 
cause  such  a  problem:  (1)  errors  in  application  of  the  boundary  conditions,  (2)  application 
of  the  boundary-layer  routine  outside  the  region  of  its  empirical  correlations,  or  (3)  the  viscous 
mixing  between  the  plume  and  free  stream  not  being  adequately  simulated  by  the  inviscid 
numerical  viscosity. 

The  boundary  conditions  used  are  those  presented  in  Fig.  8.  The  correctness  of  the 
application  of  each  was  verified.  Contour  plots  were  generated  to  ensure  transparency  at 
the  overlapping  boundaries  of  the  two  grid  sections.  The  inverse  boundary-layer  routine  in 
AR02P  includes  empirical  correlations.  The  correlations  were  not  developed  from  nozzle 
afterbody  data,  although  there  is  no  reason  to  believe  that  they  are  not  applicable.  The  effort 
required  to  investigate  and/or  expand  the  correlations  is  beyond  the  scope  of  the  present  effort. 
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Figure  14.  Variation  of  pressure  distribution  with  number  of  time  steps  per  boundary-layer 
call  (NV). 

There  are  many  ways  to  approach  the  plume  mixing  problem.  For  this  investigation,  a 
major  objective  was  to  keep  the  correction  logistically  simple,  i.e.,  avoid  methods  that  require 
the  mesh  to  change  with  free-stream  or  jet  conditions,  and  be  applicable  to  a  wide  range 
of  nozzle  afterbody  problems  with  little  modification.  The  approach  dictated  that  the  correction 
be  applied  along  the  boundary  where  the  two  grid  sections  meet  (Fig.  8).  The  obvious  parameter 
to  use  as  a  correction  for  plume  mixing  in  AR02P  is  gv  (density  x  velocity,  normalized 
by  free-stream  density),  which  is  the  parameter  used  to  correct  for  the  viscous  effects  along 
the  body.  Initially,  a  uniform  distribution  of  gv  along  the  plume  boundary  was  used.  To 
determine  what  values  of  gv  would  be  required  to  predict  the  correct  variation  of  drag  with 
NPR,  solutions  were  calculated  to  produce  pressure  distributions  which  yielded  the  same 
integrated  drag  values  as  available  experimental  data.  Figure  16  presents  the  calculated  and 
measured  pressure  distributions  at  Mach  number  0.6  and  design  NPR  of  4.1 .  Albeit  probably 
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Figure  IS.  Experimental  and  computational  variation  of  NPR  with  afterbody  drag  for  the 
inviscld  and  viscid/invisdd  Euler  code. 
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figure  16.  Experimental  versus  computational  pressure  distribution  at  design  NPR. 


fortuitous,  this  solution  requires  no  correction.  Figure  17  shows  the  results  at  an  NPR  of 
10.0.  The  solid  line  shows  the  calculated  CP  distribution  with  no  mixing  correction.  The 
dflsiipH  line  was  calculated  utilizing  a  qv  value  of  - 1 .4,  which  provides  a  satisfactory  match 
with  the  data.  Several  Mach  number  0.6  conditions  were  calculated  for  different  NPR’s  in 
the  same  manner.  The  values  of  qv  required  to  match  the  experimental  drag  data  are  shown 
as  discrete  data  points  in  Fig.  18.  The  solid  line  represents  a  modification  of  the  prediction 
technique  for  jet  entrainment  developed  by  Bauer  {Ref.  12).  Bauer’s  technique  uses  free- 
stream  and  jet  exhaust  conditions,  including  the  maximum  inviscid  plume  diameter,  along 
with  minimal  geometry  information  to  provide  a  coarse  theoretical  estimate  of  the  induced 
velocity  in  the  mixing  region.  The  induced  velocity  is  then  used  to  make  a  first-order  estimate 
of  the  effect  of  jet  mixing  on  afterbody  drag.  The  modification  was  produced  by  converting 
Bauer’s  method  to  a  density  x  velocity  term  such  as  used  in  AR02P.  The  increment  of  qv 
values,  taken  over  the  desired  NPR  range,  was  then  scaled  to  match  the  range  of  the  AR02P 
solutions  and  produced  the  rather  good  correlation  shown  in  Fig.  18.  Utilizing  the  modified 
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Figure  17.  Variation  of  pressure  distribution  with  plume  mixing. 


Bauer  method,  a  prediction  was  made  of  the  required  qv  values  to  simulate  jet  entrainment 
at  Mach  number  0.9  shown  by  the  upper  solid  line  in  Fig.  18.  AR02P  calculations  made 
to  evaluate  the  predictions  for  Mach  0.9  are  presented  in  Figs.  19  and  20.  The  result  of  using 
the  predicted  pv  value  from  the  modified  Bauer  theory  is  shown  as  the  dashed  line  on  each 
plot;  obviously,  the  theory  is  unsatisfactory.  The  solid  lines  were  obtained  by  increasing  the 
qv  terms  to  obtain  a  more  reasonable  match  to  the  experimental  data.  Even  with  the  higher 
qv  values  the  comparison  is  only  marginally  satisfactory,  although  the  integrated  drag  values 
are  the  same.  Figure  21  shows  the  values  of  qv  required  in  AR02P  to  match  the  Mach  number 
0.9  drag.  The  Bauer  theory  predicts  smaller  values  of  qv  for  Mach  number  0.9  than  for  Mach 
number  0.6,  when  in  reality  much  larger  values  of  qv  are  required  to  produce  marginal 
agreement  with  experiment. 

Past  performance  of  the  inviscid  routine  from  AR02P  has  shown  that  given  the  right 
boundary-layer  correction,  it  can  calculate  the  correct  pressure  distribution,  as  presented  in 
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Figure  19.  Variation  of  pressure  distribution  with  plume  mixing. 

Fig.  7.  Therefore,  AR02P’s  inability  to  compute  satisfactory  pressure  distributions  on  the 
AGARD  15-deg  boattail  afterbody  at  Mach  number  0.9,  even  with  the  large  values  of  @v, 
indicates  the  boundary-layer  simulation  is  not  correct  for  this  application;  however,  it  is  not 
apparent  if  the  problem  is  caused  by  applying  the  boundary-layer  routine  outside  its  correlation 
range  or  if  the  mixing  region  between  the  plume  and  free  stream  is  not  modeled  properly. 
The  manner  in  which  the  code  misses  the  minimum  pressures  as  shown  in  Figs.  19  and  20 
is  the  result  of  a  too  thick  boundary  layer  in  that  region.  Figure  22  shows  the  result  of  arbitrarily 
halving  the  boundary-layer  correction.  The  slope  of  the  recompression  region  looks  better, 
as  does  the  minimum  pressure,  but  overall  the  comparison  is  worse.  If  the  difficulty  is  in 
the  boundary-layer  routine,  it  is  not  uniform. 

The  possibility  also  exists  that  the  input  information  provided  to  the  boundary-layer  routine 
by  the  inviscid  portion  of  the  code  is  incorrect  because  of  poor  simulation  of  the  mixing 
process.  Calculations  were  made  to  investigate  the  effect  of  the  shape  of  the  gv  distribution 
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Figure  20.  Variation  of  pressure  distribution  with  plume  mixing. 


being  applied  along  the  plume  grid  boundary.  Figure  23  presents  the  variation  in  pressure 
distribution  for  two  modifications  of  the  gv  distribution.  The  tapered  distributions  were  applied 
along  the  plume  grid  boundary  from  the  trailing  edge  to  a  point  downstream  of  the  maximum 
plume  diameter,  as  calculated  by  the  method  of  characteristics.  Applying  the  modified  gv 
distributions  does  not  produce  the  desired  pressure  distribution. 

Although  the  code  calculates  excellent  afterbody  pressure  distributions  with  a  solid  plume 
simulator,  as  presented  in  Fig.  7,  the  addition  of  the  plume  and  boundary  layer  to  the 
calculations  did  not  produce  the  desired  results.  The  calculated  variation  of  drag  with  NPR 
is  in  the  right  direction  but  the  magnitudes  are  not  satisfactory.  Attempts  to  model  the  jet 
mixing  in  the  inviscid  code  with  a  simple  theory  failed  to  adequately  correct  the  problem 
so  that  the  code  could  make  useful  predictions.  No  simple  modification  was  found  to  increase 
accuracy. 
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Figure  21.  Mixing  correlation  prediction  resnits. 
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Figure  22.  Pressure  distributions  for  uniform  variation  in  inverse  boundary  layer  qv  predictions. 


3.0  CONCLUDING  REMARKS 

A  CFD  code  was  developed  which  coupled  an  inverse  boundary-layer  routine  with  an 
Euler  code  incorporating  a  jet  plume.  Calculations  were  made  on  the  axisymmetric  AGARD 
15-deg  boattail  afterbody  for  variations  in  NPR  at  Mach  numbers  0.6  and  0.9,  and  compared 
to  experimental  results.  Although  the  code  predicted  drag  changes  in  the  right  direction,  the 
accuracy  was  poor.  An  attempt  was  made  to  incorporate  a  relatively  simple  model  of  jet 
mixing,  but  the  results  indicated  that  a  more  in-depth  understanding  of  the  jet  mixing  process 
is  required  to  produce  such  a  modification. 

To  satisfy  the  needs  for  a  three-dimensional  nozzle  afterbody  analysis  tool  within  the 
current  state  of  the  art,  it  appears  that  a  Navier-Stokes  code  is  required.  Jacocks,  et  al.  (Ref. 
4),  demonstrated  the  capabilities  of  such  an  approach  for  axisymmetric  flows.  Successful 
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Figure  23.  Variation  in  pressure  distribution  with  mixing  parameter  distribution. 


application  of  a  Navier-Stokes  code  for  a  three-dimensional  situation  would  be  very  useful 
in  providing  insights  into  the  physics  of  the  complex  boundary-layer/plume  interaction 
phenomena.  Perhaps  one  can  learn  from  such  a  code  an  acceptable  means  of  representing 
the  afterbody  boundary  layer  and  plume  mixing,  thereby  allowing  the  use  of  a  viscid/inviscid 
approach. 
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NOMENCLATURE 

A/A*  Nozzle  exit  area-to-throat  area  ratio 

Cd  Afterbody  drag  coefficient  based  on  pressure  integrations 

CFL  Courant-Friedrich-Lewy  number 

CP  Afterbody  surface  pressure  coefficient 

Mw  Free-stream  Mach  number 

NPR  Nozzle  total  pressure  to  free-stream  static  pressure  ratio 

NV  The  number  of  inviscid  iterations  between  boundary-layer  corrections 

OMEGA  or  12  Smoothing  parameter 
R  Body  radius 

Re  Free-stream  Reynolds  number,  per  foot 

v  Velocity  in  y-direction 

x,y,z  Cartesian  coordinates 

X/L  Nondimensionalized  body  length 

0n  Nozzle  divergence  half-angle,  deg 

e  Density 
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